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Abstract 

We have calculated the 0(aa s ) contributions to the relationship between the MS- 
mass and the pole of the t-quark propagator in the Standard Model in the limit of 
a diagonal CKM matrix and for a massless b-quark. Analytical results for the so far 
unknown master-integrals appearing in the calculation are also given. 

1 Introduction 

Electroweak precision observables play an important role in the verification of the Stan- 
dard Model (SM) at the quantum level. By comparing precision data with corresponding 
predictions it is possible to get constraints on unknown parameters, like the Higgs boson 
mass, or unveil new physics. For electroweak processes perturbative calculations work and 
converge rather well. However, as the complexity of such calculations grows dramatically 
with the order of the perturbation expansion, complete higher order results are available in 
a few cases only. Therefore, perturbative calculations are possible only at limited accuracy. 
We usually distinguish two types of uncertainties: unknown higher order effects and the 
experimental 1 errors of the input parameters. For the calculation of electroweak observables 
the generally accepted renormalization scheme, denning a particular parametrization, is the 
so called on-shell scheme pQ [H], where, in addition to the fine structure constant (and/or 
the Fermi constant), the pole masses of particles serve as input parameters. Quark masses 
require special consideration in this context, since on-shell quark masses are not accessible 
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experimentally. Fortunately, the light quark masses in high energy processes often can be 
neglected (effects 0(m 2 /M^)) and thus can be treated as massless in practical calculations. 
The top quark is different. The large numerical value of the top quark mass in conjunction 
with the violation of the Appelquist-Carazzone theorem [H] as a consequence of the Higgs 
mechanism of mass generation, implies that a class of radiative corrections are proportional 
to positive powers of the top-quark mass which gives rise to sizeable effects. Moreover, 
the concept of a pole mass of a quark is intrinsically ambiguous due to strong interaction 
renormalon contributions [7j. 

The present experimental error of the top-quark mass, ~ 5 GeV, will be reduced to 1-2 
GeV at the LHC and/or at a future Linear Collider. The question about which type of 
mass will be determined actually with this accuracy is discussed in jS]. The higher order 
theoretical uncertainties may be estimated from scheme- and scale-variations [9 . 

The aim of the present paper is the analytical calculation of the 0(aa s ) correction to the 
relationship between the pole- and the MS-mass of the top-quark. This also provides us the 
two-loop on-shell mass counter-term for the top quark. The one-loop results of order 0(a s ) 
and 0(a) have been presented in [TU] and e.g. in jl] 2 , respectively. The one-loop 0(a) 
correction to the relationship between the top- Yukawa coupling and the pole mass of the 
top-quark has been calculated in [TT]. The two-loop 0(a 2 s ) correction is given in JIZ] . and the 
same order result obtained via regularization by dimensional reduction may be found in ^3] • 
The renormalized off-shell fermion propagator of order 0(a 2 ) has been worked out in [T3] . 
Only recently, in [15] . the three-loop 0(a^) correction has been published. Finally, the two- 
loop 0(aa s ) and 0(a 2 ) corrections have been calculated in the approximation of vanishing 
electroweak gauge couplings [TH] . Our calculation, presented here, extends previous two-loop 
0(aa s ) calculations of the gauge boson self-energies [T7] and the SM 0(a 2 ) corrections to 
the relation between the pole and the MS masses of the gauge bosons Z and W, presented 
in [HI CH]. 

The paper is organized as follows: in Sec. 2 we outline the calculation of the on-shell 
fermion self-energy. The relevant master-integrals are presented in Sec. 3. In Sec. 4 we 
discuss the renormalization problems at two-loops. Sec. 5 contains the main results of our 
calculation and in Sec. 6 we summarize our conclusions. 



2 The quark pole mass: definition and calculation 



The definition of the top-quark pole mass has been discussed in [2D] . Starting point of our 
consideration is the tensor decomposition of the one-particle irreducible self-energy of a 
massive fermion S(p, m, . . .) which, within the SM, has the form 



E(p, m, 



ip 



A(p 2 ,m, . . .)-7 5 C(p 2 ,m, . . .) 



+ m 



B(p 2 ,m,...)- l5 D(p 2 ,m,...) (2.1 



where A,B,C,D are Lorentz scalar functions depending on all parameters of the SM. For 
simplicity of the notation we will indicate explicitly only the external momentum p 2 and the 



2 (see also Eq. (B.5) in Appendix B of [TT)]') 
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mass m of the quark under consideration. The 0(aa s ) correction has a structure similar to 
the one in QCD, where C = D = 0. In this case, the position of the pole — M is defined as 
the formal solution for ip (in Euclidean metric [21] with p = jp, p 2 = p 2 \ on-shell: p 



p = — to ) at which the inverse of the connected full propagator equals zero. Thus 



ip + m — S(p, to, . . .) = ip (l — A(p 2 



+ TO 



(l - B(p 2 







im, 



(2.2) 



for p = iM. The immediate question arising here is, what is the interpretation of the complex 
mass M = M' — ~ V ? In general in a perturbative calculation a transition amplitude T 
exhibiting an unstable particle resonance has the form (U, V some spinor valued amplitudes) 



U 



1 



V 



U {-ip + M) V 



R 



ip + M (l + 0(p 2 + M 2 )) p 2 + M 2 + 0((p 2 + M 2 ) 2 ) p 2 + M 2 



+ B{p 2 ) (2.3) 



where R is a complex number and B is a background complex scalar amplitude which is 
regular at p 2 = —M 2 . T thus has the same general form as in the case of a bosonic resonance. 
We thus define the pole mass M and the on-shell width T as in the bosonic case by 



M 



M l - iMT = M' 2 - T" 2 /A - iM'V 



V2 



such that 



M 



M' 



M' 2 -T' 2 /A ■ r = — r' 



(2.4) 



(2.5) 



Since M = M' + 0(a 2 ) and r = T' + 0(a 2 ) for the 0(aa s ) terms considered in this paper 
we can identify M = M' and r = T' in the following. 

For the remainder of the paper we will adopt the following notation: capital M ~ Re M 
always denotes the pole mass; lower case to stands for the renormalized mass in the MS 
scheme, while too denotes the bare mass. The on-shell width is given by T ~ —2 Im M. In 
addition we use e, g and g s to denote the U(l) em , S77(2)l and SU(3) C couplings of the SM 
in the MS scheme. 

In perturbation theory ()2.2j) is to be solved order by order. For this aim we expand the 
self-energy function about the lowest order solution p = im§. 



E(p, m, . . .) = S \ p=imo + (ip + mo) 
and define dimensionless "on-shell" amplitudes S, S' by 



+ 



p=imo 



and 



p=imo 



p=imo 



\d(ip)J 



—moA + toq-B 



= -m S(TO , . . .) 



(2.6) 



(2.7) 



p=imo 



A + 2p 2 A + 2m\B 



= S'(to , ■ ■ .) 



(2.8) 
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where X(p 2 , . . .) denotes the derivative of X(p 2 , . . .) with respect to p 2 . To two loops we 
then have the solution (in agreement with [2*U] ) 

M , , , 

— = l + Sx + Sa + SxSi 2.9 
m 

where is the bare (m = mo) or MS-renormalized (m the MS-mass) L-loop contribution 
to the amplitudes defined in ()2.7|) and ()2.8|) . 






Figure 1: The two-loop one-particle irreducible diagrams contributing to the pole mass of 
a quark. cj) and are the neutral and the charged pseudo-Goldstone bosons, respectively. 
The number of diagrams is 24. 

According to Eq. ()2.9|) we need to calculate propagator-type diagrams up to two loops 
on-shell. The set of two-loop one-particle irreducible diagrams is shown in Fig. ^ In order 
to get manifestly gauge invariant results the Higgs tadpole diagrams, shown in Fig. |2] on 
the left side, should be included As was demonstrated in ^H] for the self-consistency 
of the renormalization group (RG) approach the two-loop tadpoles of the type shown on 
the right side of Fig. |21 should be included as well. For each quark species, there is one 
0(aa s ) two-loop tadpole diagram which is gauge invariant. For a quark with mass m q its 
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Figure 2: The two-loop tadpole type diagrams which should be included for manifest gauge 
and renormalization group invariance. 



bare contribution to the location of the pole of the top-quark propagator reads (see Sec. 4.3 
in HU) 

A t = -J-J^ im ^ 1 {d-l) r4 

167r 2 167r 2 ' (d - 4) 2 (d - 3) m w K q! 

with Cf = 4/3 in the SM and N c is the number of colors (N c = 3). The one-loop 0(a) 
result for the pole mass of a quark in the SM, in the approximation of a diagonal Cabbibo- 
Kobayashi-Maskawa matrix is given by Eq. (B.5) in Appendix B of PS]. For the calculation 
of the two-loop propagator type diagrams with several mass scales we will use Tarasov's 
recurrence relations [22] which allow us to reduce all diagrams to a few master-integrals. 
The package ONSHELL2 [23] is used for the calculation of the single scale diagrams. 



3 Master— integrals 

This section is devoted to the calculation of the so far unknown two-loop master-integrals 
needed for our calculation and shown in Fig. El We denote all master-integrals by Tab---, 
where the first letter T = F,V, J indicates the topology in accordance with the notation intro- 
duced in (22j; indices A, B, ■ ■ ■ — 0, 1, 2 characterize the relation between the corresponding 
internal mass to the external momentum: indicates a massless line, 1 corresponds to "inter- 
nal mass equal to external momentum" and 2 means that mass and momentum are different 
(see Fig. El for details). In our normalization each loop is divided by (47r) 2_e r(l + e). We 
will also use the short notations 

n- n r d n k 1 d n k 2 

Jmi ™ " p ( 3 _ « \j [{kl _ pj2 + ^ [k l + m 2] [{kl _ h) 2 + m 3] 
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7T 



-n/2 



d n k x 



AM 



n-2 



r ( 3 - s i k\ + M 2 (n-2)(n-4) 



^o(M) 



for the auxiliary integrals appearing in our calculation. 



M 

£3 




(3.10) 
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Figure 3: The master diagrams arising in this two-loop calculation. Bold, thin and dashed 
lines correspond to off-shell massive, on-shell massive and to massless propagators, respec- 
tively. 



3.1 J, 



012 



The analytical result for this type of integral for arbitrary values of the internal masses, 
momentum and powers of the propagators can be calculated by using the Mellin-Barnes 
technique j2S]- It is expressible in terms of Appell hypergeometric function F 4 [26j: 



Joi 2 (cx,/5,a,p 2 ,m 2 ,M 
= (M 2 )™~ CT ~ a-/3 



d n k 1 d n k 2 



r(f 



((p - h) 2 ) a {{h - k 2 f + M 2 ) a {k 2 + m 2 Y 
o) { 



r(a)r(a)r(/3)r(f)p( 3 -f)l 

r (§-/?) r {a+f3+a-n) V (/3+a-%) F 4 (a+(3-^,a+(3+a-n; f , l+/3-f 



2 2 

— p m 



+ l^j " " T I»r (a+a-f ) F 4 (a+a-f , a; f , 1 + 



(n/2-/3) 



-p 2 m 



M 2 ' M 2 

2 \ "1 



M 2 'M 2 



>(3.11) 



The result is symmetric with respect to the exchange M 2 <-> m 2 , a <->• (3. The e-expansion 
for this diagram up to the finite part has been given in and was recently recalculated 
in j2H]- In our case we need also the term linear in e, however, only for the particular case, 
when the external momentum is on-shell with respect to one of the internal masses. In this 
case, taking p 2 = —m 2 (on-shell index /3), we simplify our notation and write 



Joi2(v,l3,a) = J i2((T,P,a,p 2 ,m 2 ,M 2 ) 
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We thus consider 

J 012 (a,P,a) = (M 2 r-°- a ^ 



r( ff )r(a)rC9)r(f)r»(3 - f ) 

r (2-/3) r (a+/3+a-n) T (/3+cr-f ) 4 F E 



a+P+a-n, /3+a-f, §, ±±2 



+ 



m 
M 2 



r(/3-§)r( < T)r(a+ < 7-2) 4 F 3 



(X^ I 0" 2 5 - ) * ■■> 
1 + 2-/3,71-/3 2 



Am 2 

Am 2 \ 
M 2 " J 



(3.12) 

2 y 1 ' v r, 2 ivi" j j 

in the following. Let us remind that for the class of integrals J012 there are three master- 
integrals of this type: Joi2(l, 1, 1), </oi2(l, 2, 1) and Joi2(l, 1, 2) [22]- However, other indepen- 
dent combinations happen to be more suitable for performing the e-expansion: Joi2(l, 2, 2), 
J 012 (2, 2, 1), and [J i 2 (l, 2, 2) + J 012 (2, 1, 2) + J 012 (2, 2, 1)] (see also in Ref. [HI). The latter 
combination corresponds to the integral Joi2(l> 1)1) i n 2 — 2e dimensions j3U] and we have: 



^012(1, 2, 2) 



(M 2 )- 



l-2e 



r(l - e)r(l + 2e) 



l,§,l + 2e 



e(l-e) 



r(l + e 

2 \ - £ 

3^2 



-3^2 



2,2 



4m 2 
TP 



m 



r \ 1> 1 2 ^ 
' 1 2-£,2-2£ 



4m 
M 2 " 



(3.13) 



Joia(2,2,l) 



(M 2 ) 



2\-l-2e 



(l + e)r(l-e)r(l + 2e) 



e 2 (l-e) 



m 
M 2 



r(l + e) 

2 \ - £ 



-4-F3 



l,§,l + 2e,2 + e 
2,2 -e,l + £ 



3-^2 



2,l + £,§-£ 

2 - e, 2 - 2 £ 



4m 2 
M 2 " 



4m 2 
M 2 " 



(3.14) 



J c = J 012 (l, 2, 2) + J 012 (2, 1, 2) + J 012 (2, 2, 1) = (M ) 



2\-l-2e 



r(l-e)r(2 + 2e) 



r(i + £ ) 



m 2 l+j£ / l,|,2 + 2£ 
1 + 2 M^T^7 sF2 2,2-£ 



m 
M 2 



m 2 1 + £ 

1 + 2 — 

M 2 l-£ 3 2 



1,2 + £,| -£ 

2 - £, 2 - 2£ 



1 

"2 x 

£^ 

4m 2 ' 

4m 2 \ 
M 2 / 



(3.15) 



The £-expansion for the given hypergeometric functions were worked out in ^H] and [3*T] . 
While the individual expressions are rather lengthy the results for the £-expansion of the 
integrals have the compact form 



J 012 (l,2,2) = (M ) 
+£ 



ln(l + y) 



ln(l + y) - ln(y) 
1 



2 ln d (l + y) - 2( 2 ln(l + y) - 3 In y hr (1 + y) + - hT 3/ ln(l + y) 



-2 ln(y)Li 2 (-y) - 3 ln(y)Li 2 (y) + 6Li 3 (-y) + 6Li 3 (y) 



+ 0(£ 2 



(3.16) 



•W2,2,l) 



(M 2 ) 



2\-l-2e 



(1 + 2/) 2 f 1 



(1-5) (l-y)[e 



i-y 



ln(l + y) + In y 



-hi 2 y + 2C 2 



+2- ln(l + y) ( ln(l + y) - my 



+£ 



1/ 

i + y 
y 



v ) " In y ln(l - y) + 4Li 2 (-y) + 3Li 2 (y) 

12S lj2 (y 2 ) - 4S 1;2 (-y) - 6S 1)2 (y) + 241n(l - y)Li 2 (-y) 



+ 18 ln(l - y)Li 2 (y) + 9 lny ln 2 (l - y) + - In 2 y ln(l - y) - 6(2 ln(l - y) 
+^ ln 3 (l + y) - 2( 2 ln(l + y) - 4 lny ln 2 (l + y) + In 2 y ln(l + y)) 

+^ln 3 y - 2C 3 + lnyLi 2 (y) + 4^ lnyLi 2 (-y) 

6 y y 



7^+4 . Q 9y + 5 T , 
-4 Li 3 (-y) - 3 Li 3 (y) 



y 



(3.17) 



J c = J 012 (l, 2, 2) + J 012 (2, 1, 2) + J 012 (2, 2, 1) 



(M ) 

+€ 



1 -y U 



6 In y ln(l - y) + - In 2 y - 2( 2 + 8Li 2 (-y) + 6Li 2 (y) 



24S li2 (y 2 ) -8S li2 (-y) - 12S lj2 (y) - 12( 2 ln(l - y) + 48 ln(l - y)Li 2 (-y) 

1 



+36 ln(l - y)Li 2 (y) + 18 lny ln 2 (l - y) + 3 In 2 y ln(l - y) + - In 3 y 



+201nyLi 2 (-y) + 241nyLi 2 (y) - 2( 3 - 44Li 3 (-y) - 42Li 3 (y) 



where 



y = 



M 2 

,2 



(i + y) 2 



rrr y 

The expressions for the original set of master-integrals then read: 

Jr. 



+ 0(e 2 ) , (3.18) 



(3.19) 



9 9x 9 [4m 2 -M 2 (9n-32)l 
J i 2 (l,l,l,m 2 ,M 2 ) = M 2 ' 



(3n- 10)(3n-8)(n-3) 
2(M 2 - m 2 ) [4m 2 (2n - 7) + M 2 (n - 4)] 



-A (m)A (M) 
J 012 (l, 2,1, m 2 ,M 



(3n- 10)(3n-8)(n-3) 

n - 2) 2 [(7n - 24)M 2 + m 2 (8n - 28)] 



J i 2 (2,2,l) + J 012 (l,2,2) 



4M 2 m 2 (3n - 10)(3ra - 8)(n - 3)) 
,, 2(M 2 -m 2 )(2n- 7) 



(3n-10)(n-3) 



-Joi 2 (l,2,2)- 



M 2 



(3n- 10)(n-3) 



Jr 
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\M 2 (n-A) -2m 2 (2n-7)l , s , , SA/ s (n - 2) 2 (2n - 7) 

;j -Joi 2 (2,2,l) + A (m)A (M)- 



(3n-10)(n-3) v ' ' 7 v ' v MM 2 m 2 (3n - 10)(n - 3) 



J i 2 (l,l,2,m 2 ,M 



[2m 2 (n-3) + M>-4)] 



(3n- 10)(n-3) 
(3n-ll)M 2 w (n-2) 2 (2n-7) 



J i 2 (2,2,l) + J i 2 (l,2,2) 



' (3n - 10)(n - 3) Jc + A ° WA ° (M) 4MV(3,-lQ)(,-3) ' (3 - 20) 
where A (m) is defined in (I3.1U|) . 

3.2 V0012 

Let us consider the second diagram of Fig. El 

V 00mM (K a, (3,a) = J {{p _ ki y )X{{ki _ k 2 )2)t(kl + m?Y{kl + M 2 ) a ' (3 ' 21) 

with the external momentum on the mass-shell, p 2 = —m 2 . For this class of integrals we have 
only one master-integral, the one with indices (1, 1, 1, 1). For constructing the e-expansion 
we may use the differential equation method [32J. In terms of the new variable r = m 2 /M 2 
the result may be represented as 

VoomM = (m 2 )- 2e (j- 2 + ^ + F(r) + eE{r) + 0{e 2 )^ . (3.22) 

The analytical results for arbitrary r read 

. . 5 (1 — r) ln(l - r) 

S{r) = - + - — - + lnr, 

2 r 

F(r) = y-2— — Li 2 (r)+C 2 (l-3r)+41nr+ M --Jln 2 r 

+ (1 ~ r)2 Inr ln(l - r) + 4(1 = r) ln(l - r) - (1 = r)(3 = r) ln 2 (l - r) , 
E{r) = |_(l-r)(3 r -7) Si2(r)+ (l- r) (3 r -5) L . 3(r) 

H — — -Li 2 (r) ln(l — r) Li 2 (r) In r 



2(1 -r)*, , (1-^(1-2^ , 



lnr In (1 — r) H In r ln(l — r) 

/• r 

(4-3r). 3 (l-r)(5-3r). 3/1 , „ _ 
In 3 r + i - i ln 3 (l - r) - (1 + 3r)C 3 



6 3r 
-2(1 - r) C2 lnr - In(l - r) - ^^Li 2 (r) 

2(1 - r)(r - 3) . , 4(1 - r) 2 , „,„ , , , 



ln 2 (l - r) + — — lnrln(l - r) + 2(2 - r) hTr 

/ r 

12fl — rl 

+6(1 - 2r)( 2 + — ln(l - r) + 12 lnr . (3.23) 



For r = 1 (on-shell case) we have: 

5 



< s '".) = ~, F(l) = y-2C 2 , ^(l) = f -6G-4C3- 



3.3 V1112 

The last and most complicated master-integral is 

VmmmM(o>, (3, <T, A) = 



d n hd n k 2 



(Go - £a) 2 + m 2 ) CT ((fc 1 - A; 2 ) 2 + m 2 ) a (k 2 + w?Y(k\ + M 2 ) ; 



(3.24) 



where the external momentum is on the mass-shell, p 2 = —m 2 . Let us introduce here an 
angle defined via [33 j 



M 

cos = — , M < 2m 
2m 



The e-expansion now can be written as 

V m mmM(9) = (m 2 )~ 2e 



^ 2 + - £ S{9) + F{9)+eE{9) + 0{e 2 ) 



(3.25) 



(3.26) 



The analytical results for arbitrary values of the angle read 
5 



S{9) 
F{9) 

E{9) 



- AO cos 6 sin 9 - 4 cos 2 9 In (2 cos 9) , 
19 

— + 49 sin(20) In (2 sin 20) - 89 sin(20) - 40 2 sin 2 9 + 2( 2 sin 2 9 



sin(20) 



4Ls 2 (20) - 3Ls 2 (49) 



+ 4 cos 2 9 In 2 (2 cos 0) - 16 cos 2 9 In (2 cos ( 



65 



- 249 sin(20) - 160 2 sin 2 9 + 12C1 3 (20) sin 2 9 - 3C1 3 (40) sin 2 9 



+169 sin(20) In (2 sin(20)) + 160 2 sin 2 In (2 sin 0) - 40 sin(20) In 2 (2 sin(20)) 
+8 sin(20) In (2 sin(20)) Ls 2 (20) - 6 sin(20) In (2 sin(20)) Ls 2 (40) 
-48 In (2 cos 0) + 80 2 sin 2 In (2 cos 0) + 48 sin 2 In (2 cos 0) 

+ 16 cos 2 In 2 (2 cos 0) - - cos 2 In 3 (2 cos0) + 80 sin 2 0Ls 2 (20) 

3 

-16sin(20)Ls 2 (20) + 12sin(20)Ls 2 (40) - 4sin(20)Ls 3 (20) + 3sin(20)Ls 3 (40) 
-12Ca sin 2 In (2 cos 0) - 4( 2 + 8( 2 sin 2 - 2( 3 sin 2 , (3.27) 



where the Lsj (0) are so-called log-sine integrals jHS] defined by 



J d<p ln^ 1 


2 am- 
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(3.28) 



Clj (9) are Clausen functions, 



CI,- (9) 



2i 

1 
2 



Hj (e ie ) -Lij (e" ie ) , j even 
Li, (e^+Lij- (e~ ie )l , j odd 



(3.29) 



and the Lij are poly-logarithms. 

For M 2 = m 2 the angle 9 is equal to 7r/3 and the integral is equal to the single scale 
diagram Villi calculated in [36 . In this case we obtain 



3 



7T 



E 



2 v/3 
65 



7T 



19 7T 



7T 



--^=, F - = — + -=ln3-4-=--C 2 -7 



7T 



6C 2 + 4^1n3 



1 7T 



9 



V3 
Ls 2 (f ) 



2^/3 



In 2 3-12 



V3 2 

7T 9 7T 

V3 ~ 2^3 



"2 C3 + 7 ^3 



In 3 - 2! 



v/3 



+ -7rLs 2 ( - 



L *» (f ) 

V3 



C 2 + 4C 2 ln3 

2lMf) 
2 73 



(3.30) 



Another case of interest is the massless case M 2 = 0, which corresponds to 9 = ir/2. In 
this case the diagram can be reduced to J mmm (l, 1, 1) jHZllSH! plus a combination of simpler 
vacuum diagrams JHHI- The results here are 



7T 



7T 



19 



7T 



65 



The expression ()3.27|) is directly applicable in the region M < 2m only. For (M > 2m) one 
has to perform the proper analytical continuation. How this can be done is described in 
details in |HJ |MJ EI] • ^ or ^ ms P ur P ose ; let us introduce the new variable 



y = e 



Jo-26 



ln(— y — iaO) = \ny — ian. 



(3.3i; 



which coincides with the variable y defined in f)3. 19|) . In terms of this variable the expression 
(I3.27|) may be written as 



S(y) 
F(y) 



5 (1 4- v) 2 

^ + (l + y)lny-^-^-ln(l + y) 
2 y 



19 l-y 2 
— + — 

2 y 



Li 2 (j/) + 3Li 2 (-y) 



+ ln 2 y+(l-y)C 2 



Ulny - 21ny ln(l - y) - 21nyln(l + y) 



y 



hr(l + y) - 41n(l + y) + lnyln(l - y) 



J5(y) = tt 



65 l-y 2 
- + — 

2 y 



4S li2 (j/) - 3Si, 2 (y 2 ) - 6 ln(l - y)Li 2 (-y) - 2 ln(l - y)Li 2 (y) 



11 



+4Li 3 (y) - 6 ln(l + y)Li 2 (-y) - 2 ln(l + y)Li 2 (y) - 2 lny ln(l - y) ln(l + y) 
- Iny ln 2 (l - y) - C2 ln(l - y) + 12Li 2 (-y) + 4Li 2 (y) + 4 lny ln(l - y) 



-12C 2 ln(l + y) +41n 2 y + 2. 



l + 2y)(l + y) 



C 2 hi(l+y)-41n 2 yln(l-y) 



-2 In 2 y ln(l + y) + In 3 y + (1 - y) (2( 2 lny - Cs) + 4y In yLi 2 (y) - 4y( 2 



+(1+2/) 

, (l + y) 



12 lny - 81nyln(l + y) + 21nyln 2 (l +y) ln 3 y 

3 



4 ln 2 (l + y) - 12 ln(l + y) - - ln 3 (l + y) + In 2 y ln(l - y) 



-lnyLi 2 (y) 



--(y-l)(y + 3)Li 3 (-y) . 

y 



(3.32) 



The results (|3.16j) - (|3.2(J|) . (|3.23|) (|3.32|) have been checked by a heavy mass expansion |2I| 
with the help of the packages described in |4D] . 

4 Renormalization 

The mass renormalization constant Z t in the MS scheme at two loops may be written in the 
form 

m ,, = mt (S) z, = mt{ /) (i + ^ j^'" + ^ ^S'" 

+ ^cr is? + + 0(9 *' a - 2) ) • (433) 

where a s = y 2 /47r. Note that in contrast to the corresponding definition for bosons (used, 
e.g., in [TSJ EH]), the mass renormalization constant here (i.e., for fermions) relates the 
masses and not the squared ones. The coefficient Z^ ,l > may be extracted from Eq.(4.31) 
of [1^] (2Z2'^ = ^iquark)- The coefficient Z^f^ is well known 10J. For completeness we 
give them here: 



1 1 m| 



rrir 



3m 2 H 



3m 2 w 



Am 2 w m 2 H 



m 



w 



3m^ 
2 m 2 H 



+ 



3 m. 



+ N C 



mi 



m 



w 



i i > 



3 

-3C/. (4.34) 
In our calculation we obtained the two-loop renormalization constants Zjffl and Z^ 2 ^ 



m 



2 9 m% mf 9 m 2 „ 9 m 2 9 m 2 

— 9./V, 



2N r . 



+ 



4 m? A/ m 2 



+ 



H 



m. 



b W" u H 

3 m, 2 
+ ^^ + 



t 

m 2 H m l w ' 8m w 



+ 



IV 



2 m| 



- 1 



19 m| 



my^mjj 2 mw 48 m 



31 

+ 777 



k ' 24 



,(4.35) 
(4.36) 
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where, in the SM, Cf = 4/3, N c = 3 and the first five quarks are treated as massless. 
The terms proportional to m\ are coming from the tadpole contribution and will cancel in 
observable quantities. We may use the SM renormalization group equations to cross-check 
the 1/e 2 - and l/e-terms (see also [HHHI]). The coefficient ()4.35|) of the 1/e 2 term may be 
calculated from the relations 



8 

2 u 2( 2\ 



2 * 







dmj 



1 V« d z {2) 



(4.37) 



where we adopted the notation m t n = m t 9 j = 9,9s- The anomalous 
dimension of the top-quark mass jt is defined by 



_ 1 2 & I 2\ 

7t = W mt{[l ] 



_d_ 

d 9j 



(4.38) 



Translated into a relation for the coefficients defined in (I4.33|) we have 



It 



16tt 2 



9 2 Zt 1) +9 2 s Z^) + 



2„2 



9 9 



d 



dm 2 



j (16tt 2 ) 2 
z (i,U = Z (i,i) ( z (i,D 



2^ + 0(9 



4 „2\ 



3 m? 

_ * +2 N C 



m 



m W m H, 



(4.39) 



Note that ()4.39|) explicitly reveals, that the systematic inclusion of the tadpole contributions 
is important for the self-consistency of the RG equations. 

The terms proportional to 1/e may be deduced from the RG equations calculated in the 
unbroken phase. It has been shown in fSJ (details are given in 021) that in the MS 
scheme we may write 

m ^ 2) = l^WT mV) ' (440) 

where m 2 and A are the parameters of the symmetric scalar potential and Y t is the top- 
quark Yukawa coupling. As a consequence we get the following relation for the anomalous 
dimension of the mass of the top-quark j t 



1 lPx 



(4.41) 



where the relevant RG results may be found in |4*3j : 

1 



1 2 9 2 

m 



m? d(j? 



16vr 2 



X + 3Y t -~g --g 



y2 2 



(16tt 



7y = — [i 



d 



Yt dfi 2 



Y 



16vr 2 



9 y2 _ 4 2_ 9 2 _ 17 /2 
4 * 9 S g9 



(16tt 2 ) 2 
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_d_ 



-A 



1 



(16tt 2 

+ 7 



2A 2 + Q\Y? 



9 2 s Y t 2 



167T 



2\2 



40A - 96Y t 2 



2 ° 2 
i 



Finally, the parameter relations 



Y? 



2mi 



A 



3m 2 H 



Am 2 w 



,'2 



4(m| - m 2 ^) 



(4.42) 



provide the necessary bridge between Eqs. (|4.41|) and ()4.42|) and our Eqs. ()4.34|) and ()4.36|) . 

We now turn to the discussion of the pole-mass relation ()2.9|) . The calculation of the 
one-loop MS renormalized on-shell amplitude Si (see ()2.7j) for the definition) is simple. We 
get it by rewriting the bare expression in terms of MS parameters. In terms of the amplitudes 
X defined by 



Si, 



2 2 

9s,o v (i) | 9o 



16n 



2 <Xafl 



16tt 



2 a,0 



(4.43) 



we obtain 



Si 



MS 



lim 



1 i 9 *fl y(l) , go y(l) 



167T 



167T 2 " 



o,0 



lim ( — 



*L X« + - 

4tt Qs 16tt 2 sin 2 6 W 



x« + v -z^ 



+ 



167r 2 sin 6 W 



I 1 7 (1,1) 



X (D 



(4.44) 



where X« = ACf+Zgj 1 ) In ^ and X« = ixg where is given in Eq. (B.5) of Ref. [Hj]. 
We note that the explicit /i-dependence is given by the following structure: 

,2 



X. 



(1) 



In 



mi 



(4.45) 



where AJQ does not explicitly depend on /x. At the two-loop level, we may avoid the 
consideration of the wave-function renormalization as well as the renormalization of the ghost 
sector and of the gauge parameters if we look directly at the full two-loop MS renormalized 
on-shell amplitude. The latter can be written in the form 



Ea + SiSil = lim(s 2;0 + Si, S / 10 + ^ 6 2 

d 



1^(2,1) 1 J_^(2,2) 
e aa s ' £ 2 aa s 



a: 



1 



47r 167T 2 sin 2 W £ 



2m 



1 dmlo 



y(i) , 7(1,1) 



1 + 2m 2 



<9 



y(l) 



a* 



47rl67r 2 sin 2 ^V aas /i 2 



C (2,2) ln 2^ + C <(2,l ) ln ^? 



14 



1 



+-ln 1 

8 V u t 



(18u; 2 + 21u t + 17) 



+ \nu t 



22u t + 17 



(1 + 2^)^/ 



1\ , 1 + £J t - £J? 2 



^ ■ ( 1+ ^^ + 7.,-9)Li 2 (I) 



4u t 



-^(1 W(l + 2u*){|s 1>2 ( 



— ) - ^Li 3 - lnu t U 2 (— 

u t J 2 \u t J \u t 



In 2 u t + 6( 2 



(l + yz) 2 (l + yl)(17 + 41yz + 17y|) 



3 2Li 3 (yz) + Li 3 (-yz) - 3( 2 ln(l + y z ) 



-2 lnyz (2Li 2 (y z ) + Li 2 (-y z )j - In 2 y z (ki(l - y z ) + ^ ln(l + y^ j 



| (l-yz)(l+yz) 3 (17 + 41yz + 17y|) x 



2 In yz ^ln(l - y z ) + i ln(l + j/ z )^ + 2Li 2 (y z ) + Li 2 (-yz) | 



4(l + yz) / 5 _ 4 wrf/z 



9 y 2 \^ + vz) 2 

(l + yf)(l + 4yz + yf) 



x 



(i + yz) 



3 2Li 3 {y z ) + Li 3 - 3( 2 ln(l + y z ) 



-2 lnyz \2U 2 (y z ) + Li 2 (-y z )j ~ In 2 yz (Ml - yz) + - ln(l + yz) 



+(l-yz)(l + 4y z + y|) 



2 lnyz (ln(l -yz) + \ ln(l + y z ) J + 2Li 2 (y z ) 



9 \ , (l + 2yz-24y| + 2y| + y|) T , 

- yz) (1 + yz) ln(l + yz) + ^ _ ^ Li 2 (-yz) 

3 yz(l + 3yz)(4-yz + yl) , 1 yz(2 + 9yz + 3y| + 16y| + 6y|) 2 
+ 4 (T^z) ^'"i (l + yz)d-yz) M ^, 



(l + yz) 3 (9 + 32yz + 9y|) T 447 125 (1 + y 2 ) 
5 Ll 2 i 2/zJ + -77- H 



32 



4^(1 - y z )y 2 z 
Vz 



1 - u t 



C 3 -4C 2 ln2 + 



16 9 yz 
(l + yz) 4 (17y|-19yz+17) 
8^y| 



ln(l +yz) 



15 



— - In 2 y z < 



685 17 67 

+ 7T7T^ + 



36 36yf 2Ay z 



335 497 2 17 3 
~2A Vz ~ ~72 Vz ~ V2 Vz 



25 



(1 + 3/z) 



51y| + 113y| + 134yf + 237yf + 197y z + 68 



_ 1 f 20 - 39j/g - yf - 80y| - 20y| 

3 2 \ yz(l-yz) 

1 [4157 425(1 + 1/1) 2561(1 + |/|) 1 
w t \ 72 72y| 96y z J 

1 . f 187 17 133 535 211 



i -yz 

lny z 



7 -49y z + \ly\- 55y| -16y|) 
(l-y z )(l + yz) 2 



, - 2 17 3 50 

+ - C ^ — + %| + 12^ + l2-^ + l2-^ + ¥^-l^ 



2/z 



(1 + yn) 1 



y H (i + yz)\ n (i + y z ) 2 



2u t (l + y H y y\ 



yz 



yu 



f iiq + yfrXi + y*) 2 , RN | i (i + yz) 4 ' 

u t {l + y H y\ Sy% +8 c+ 2 y% 



+ -C: 



9 



io t \2y H 2 4 



oo t (625 + 1286y H + 625y 2 J ) 



36 



(1 + 2/h) 2 



i (i - y H f . 

H 5 In 2/// 



y# 

112 + 82/h- lOy 2 , - 3yf, . 2 



ln(l-y H ) + -ln(l + y H ) 



(l- y 2 / )--(l + y 2 / ) hi/* 



1 1 (l + ^) 2 (5- 62^ + 5^) 



, 1 1 (1 + ^(6-63^ + 5^) , 
In y H + - my H 



Z/ff 



wio. ^ 31 nnx , (l-y g ) 2 (l + ^) 
ln(l + y H ) - 7;— C2M 1 + yff) 2 



^ (1 - y „ )(1 + y „) | (5 - 28 +5 ^) Li2 + _ 

1 (l-y g ) 2 (l + y|) [3 



UJ t 



y 2 H 



I 2 



2Li 3 (y#) + Li 3 (-y H ) 



\ny H 



2Li 2 (y H ) + Li 2 (-y#) 



(4.46) 



where 



yA 



m 2 






'-\ 


A- 




1 + \ 


A- 


4nij 



A = H,Z . 



(4.47) 



We have explicitly factorized the RG logarithms, and Cf££\ which may be calculated 

also from the one-loop result and the knowledge of the mass anomalous dimensions (see 
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for the general case): 



= Z^ 1+m* 







'dm? 



Z (l,l) _ ^(2,2) 



= 2Zi 2 i) + AZ^Z^ + ZW)AX« + Z^ I 1 + 2m; 



d 



where 







m 



i [3(i+^r 



+ - T 



8^ty 2 H 

Vh 



*<9m t 2 



(4.48) 
AX« , (4.49) 



ln(l + y H ) - 



{l + y H ){3 + y H ) 



o; t [4 y\ (l + y H y 4 



+ l(y 2 H + 3y H + l) 



Vh 



2 ^ j/ 
'1 + VhY \ 



Vh 



25 



2(1 + to) 2 1* 



+ 



1 f (1 + to) 2 19 17(1 + 1/1) 20(1+1/1) 



u t 1 3?/2 



wtj/z \ ((1 + yzY 



-5-4- 
9 V ( 

(i + yz) 3 

72u t 

169 10(l + ?/|) 1 



ln(l + yz) 



36^/1 

2/z(l + yz)(2 + yz) 



9yz 
hri/z 



,i - yz 

3 (34 + 41y z + 34y z ) ln(l + y z ) + — — — \ny z 



Vz 



72 



+ 



92/2 



+ 77 hi 1 



;i - yz) 

1 \ (l-^)(4cu 2 + c^ + 1) , 1 



+ In—. (4.50) 



The C^'s in the SM (in contrast to QCD) have non-polynomial structure in the dimen- 
sionless coupling constants which originates from the tadpole contributions. 



5 Results 

The relation between the top-propagator pole M and the MS mass m t can be written as 

— = l + {£i) +(e 2 + +0(g\a 2 s ), (5.51) 

m * 1 Jms I Jms 

where the r.h.s. is given by Eqs. (|4.44|) and ()4.46|) . 

We would like to elucidate several aspects of our calculation: all diagrams have been 
calculated in the so-called modified MS scheme (MMS) jHZ|, which is defined by multiplying 
each loop by the factor l/(47r) £ /r(l + e). It has been shown in ^Hl QEI EH] that at the 
two-loop level this scheme is equivalent to the MS scheme when applied to calculating mass 
relations. We once more emphasize that the inclusion of the tadpole diagrams, shown in 
Fig. El is important for several reasons: i) to restore gauge invariance; ii) to get UV counter- 
terms which coincide with the ones calculated in the unbroken phase; iii) to restore the 
proper form of the RG equations in the broken phase (the result (|5.51j) is manifestly RG 
invariant through 0(g 4 ,a 2 )). 
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As the top is an unstable particle the pole of the top-propagator exhibits a real part 
which is the pole-mass M t and an imaginary part which up to a sign gives the width F t 
divided by two (see ()2.5|) for the precise definitions). The imaginary part is coming from 
diagrams with W, ^ lines (see Fig. [Q) and may be calculated analytically. We find 



£i = - 2 i m j£ = * J_L 

m t m t 8tt sin Ow { 8t0i 

oc s e 2 [9, m? / , 



3 , ., 
-uj t + -LO t 



47T 87T Sin 2 6y/ 

1 



1 

4 

1 



1 / 1 \ 17 21 i£ 

8 V w < w * 



-^ln(l-c^)(l-^) 2 (— + 4) -|lnu; t (l-u t -2u, 



2+1 



out 1 2 
ln^ln(l - oj t ) + 2C 2 + 2Li 2 (u t ) 



(5.52) 



as a result in terms of MS parameters. With the help of ()5.54|) we may get the corresponding 
result in terms of on-shell parameters 




+ iln^ 

+ 2 M 2 



M 2 




M 2 



(5.53) 



which coincides with the well know result jH] for the tree and the 0(a s ) correction to the 
partial decay width for the process t — ► bW in the approximation of a massless b-quark. 

Very often the inverse of the relation (|5.5H1 is required. To that end we have to solve 
the real part of (|2.9|) iteratively for m t and to express all MS parameters in terms of on-shell 
ones. The solution to two loops reads 



Re i Sl _ ~ Re ) S2 + SlSl _ + 

I J MS I J MS 



-E( A 



in 



2^(1). 



d 



dm 2 



Rel Ex \ - (Ae 



MS 




1 2 



MSJ 









MS 



(5.54) 



i 2 =M 2 



,e=eos 



where the sum runs over all species of particles j = Z, W, H, t and 



(Am 



2^(1) 



—M 2 os 

3 16vr 2 sin 2 9 



w 



(i) 



m 2 =M 2 
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stands for the self-energy of the jth particle at p 2 = —m 2 in the MS scheme and parameters 
replaced by the on-shell ones. The values of Xj are given in Appendix B of ^Sl GUI- The 
relation ()5.54|) includes also the transition from the MS to on-shell scheme for the electric 
charge. It can be described as 

(Ae) (1) |-Re|s 1 l_= {5a bos + 5a lep + fc top + A^Lonsl^I) -^a udscb (M|)} Rel Ex 1 , 

ae I J MS I J MS 



where 



and 



<5Aa udscb (M|) = — - In- z 1 



9tt V"" /i 2 3/ 
For numerical estimations of Aa^ drons (M^) we use the results of 

A «h 5 adrons( M !) = 0.027773 ± 0.000354 ; a' 1 {M 2 z ) = 128.922 ± 0.049 (5.55) 

at M z = 91.19 GeV. 

For the 0(aa s ) contribution only the expression (|5.54J) actually simplifies to 



M t 1 7ms I Jms ^le^sin 2 ^ 



^^{^^ + X «' m *lk X «\ (5 - 56) 



where 



m 



2 —X {1) - m 2 9 AV (1) I In™* m 2 dZ "' 1} I 



mi 



dm 2 " "" /i 2 \8m^ m^rriff 



For Eq. (J5.56|) we present a semi-numerical result for /i = M t and the following input param- 
eters: a = 1/137.036, M w = 80.419 GeV, M z = 91.188 GeV and M t = 174.3 GeV, 



m t {M t 



1 ~ c f~r ~ 

7T 

3 Mf 
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47rsin 2 ^ 5 
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1 Mi 



IV 
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Mj 
Mf 
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8 M/M^r 
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IV 



1 + Vj 



InV, 
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1M| 

Im7 
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-C Qs a 

/ 4?r 4vr sin 2 6g s 
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InY, 
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; (i + r^)inr ff 



1 M 2 2 + 8Yh - 10Y| - 3Y^ 
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1 M^ (11 - 50Y H + IVY 2 ) 
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M 2 W 
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1 H 




-Li 2 (-Y H ) + (1 - Y H ) 2 Li 2 (Y H ) 



- 2U 3 (Y H )+Li 3 (-Y H ) 



-InY, 



2Li 2 (Y H )+U 2 (-Y H ) 



,(5.57) 



where Y4 is defined via the pole masses as 



Y A 



1 - 



4Mf 



1 + A 1 



4M| 



A 



H,Z 



For illustration of the numerical significance of our result we shown the two-loop 0(aa s ) 
corrections to M t /mt(mt) and mt(M t )/M t (Fig. 0]) as a function of the Higgs boson mass in 
comparison with the two- and three-loop QCD ones ^2] (for simplicity, we consider the case 
a s (/x) = a s (M z ))- In a wide range of values of Higgs boson mass (100 GeV < M H < 1000 
GeV) our result is comparable in size to the 3-loop QCD corrections. The corrections grow 
with the Higgs mass, again this is a consequence of the breakdown of the decoupling of heavy 
particles in a "spontaneously broken" gauge theory. 
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Figure 4: Electroweak 0(aa s ) correction to M t /m t (m t ) — 1 [left] and m t (M t ) / M t — 1 [right], 
in comparison with 0(a 2 s ) and 0{o? s ) QCD corrections as a function of the Higgs boson mass 

6 Conclusion 

The main result of the present investigation is the two-loop 0(aa s ) relationship between 
pole- and MS-mass for the top-quark within the SM. It is given by Eq. ()4.46|) . Numerically 
its size is comparable to the 3-loop Q CD-correction for a light Higgs boson and it increases 
for a heavy Higgs boson. As a byproduct, several new massive master- integrals have been 
calculated analytically (Sec. 3). 
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